Numerical approximation of a viscoelastic frictional contact problemApproximation numérique d'un problème viscoélastique de contact avec frottement

We consider a fully discrete scheme for a quasistatic frictional contact problem between a viscoelastic body and an obstacle. The contact is bilateral, the friction is modeled with Tresca's law and the behavior of the material is described with a viscoelastic constitutive law with long memory. We state an existence and uniqueness result for the discrete solution, followed by error estimate results. Then, we present numerical simulations in the study of a two-dimensional test example. To cite this article: Á. Rodríguez-Arós et al., C. R. Mecanique 334 (2006).     

大型结构特征值问题的Lanczos子结构并行算法

本文利用子结构和Lanczos方法,提出了大型结构固有频率与模态的并行解法。该方法在Lanczos方法的求解过程中,仅利用子结构刚度阵和质量阵并行进行凝聚,进而求得新的迭代矢量,最终求得三对角阵对应的特征值和特征向量。该算法在西安交通大学ELXSI-6400并行计算机上程序实现,计算结果表明能有效地节省计算时间和计算机的内存,为一种有效的大型工程结构动力问题的求解方法。     

Real-Space Renormalization Estimates for Two-Phase Flow in Porous Media

We present a spatial renormalization group algorithm to handle immiscibletwo-phase flow in heterogeneous porous media. We call this algorithmFRACTAM-R, where FRACTAM is an acronym for Fast Renormalization Algorithmfor Correlated Transport in Anisotropic Media, and the R stands for relativepermeability. Originally, FRACTAM was an approximate iterative process thatreplaces the L × L lattice of grid blocks, representing the reservoir,by a (L/2) × (L/2) one. In fact, FRACTAM replaces the original L× L lattice by a hierarchical (fractal) lattice, in such a way thatfinding the solution of the two-phase flow equations becomes trivial. Thistriviality translates in practice into computer efficiency. For N=L ×L grid blocks we find that the computer time necessary to calculatefractional flow F(t) and pressure P(t) as a function of time scales as N^1.7 for FRACTAM-R. This should be contrasted with thecomputational time of a conventional grid simulator N^2.3. The solution we find in this way is an accurateapproximation to the direct solution of the original problem.     

An efficient algorithm for hypersonic viscous flows

The CSCM-S algorithm proposed by Lombard et al. is a very attractive tool for solving multidimensional Euler and Navier-Stokes equations. However, it is not economical due to the use of global sweeps in the whole computational domain. In this paper we suggest a modified strategy, which combines a single-marching technique for supersonic dominated region with a multi-sweep procedure for pure subsonic and complex flowfield. The new algorithm may save significantly CPU time and is more suitable for engineering applications. The project supported by National Natural Science Foundation of China     

反演在骨生长方程参数识别中的应用

基于BFGS算法，根据自行设计的动物实验，得到在不同应力环境中，快速生长期大鼠股骨生 长与重建过程中骨密度实测数据，采用反演数值方法，获取了骨生长方程中随时间变化的 生物参数B和K. 通过正演验证，表明所建立的基于BFGS算法的参数识别方法具有较好的 稳定性和较高的识别精度，能够反演出比较切合实际的生物参数数值. 研究表明，反问 题方法是解决生命现象不可测性和未知性问题的有效手段，把反演方法应用到骨生长与重建 等生命现象的研究中，可确定、修正预设的数学模式，为数值量化骨适应生物模型的建立提 供了一条可行的途径.     

计算具有区间参数结构的固有频率的优化方法

基于区间函数的单向包含性质，把具有区间非确定参数结构的固有频率所在区间范围问题 转化成两个全局优化问题，并采用一种实数编码遗传算法求取问题的全局解. 用一种能够求 得剪切型结构和弹簧质量系统特征值范围精确解的单调分析方法进行检验. 在一 些文献中，直接采用区间数运算法则和有限元法得到结构区间刚度阵和区间质量阵，并把关 于该区间刚度阵和区间质量阵的广义区间特征值问题的特征值区间作为待求的非确定性结构 的特征值所在的区间范围，该方法易于扩大问题的解域. 算例表明，可望得到结构 固有频率区间范围的准确解.     

条形平板中裂纹识别的遗传算法

研究条形平板中平行于上下边界的单一裂纹的识别问题。给出了当一声波入射这一条形平板时，裂纹参数与反射、透射系数之间的关系式，并应用遗传算法对裂纹进行了定量识别。计算结果表明本文给出的方法具有较好的识别精度。     

部分输入未知条件下结构动力复合反演的分解算法

在充分利用部分输入已确知而部分输入未知的激励特性的基础上，提出了结构动力复合反演的分解算法，该算法从源头上消除了迭代过程中参数识别与荷载反演的相互影响，降低了问题的计算规模。对于线性参数系统，该算法不经过任何迭代计算即可一次性完成结构参数识别及荷载反演。将其与松弛法结合，可解决非线性参数系统的识别问题，与文献[4]的方法比较，其收敛速度有显著提高。     

On the solution of plane boundary-value problems of elasticity in a rectangle

An algorithm for solving plane boundary-value problems of elasticity for a rectangular domain is expounded. The algorithm is based on a complex-valued representation of the general solution to the differential equations of the plane problem and on the use of Lagrange polynomials to satisfy the boundary conditions. The algorithm can quite easily be implemented in a computer program. This is probably the simplest way of solving boundary-value problems of this class __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 1, pp. 97–102, January 2006.     

On numerical solution of dynamic problems in the theory of reinforced shells

Dynamic problems for cylindrical shells reinforced with discrete ribs are examined. A numerical algorithm based on Richardson extrapolation is developed. Specific problems are solved, and the results are analyzed __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 5, pp. 50–56, May 2006.